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The electric charge density in mesoscopic superconductors with circular symmetry, i.e. disks 
and cylinders, is studied within the phenomenological Ginzburg-Landau approach. We found that 
even in the Meissner state there is a charge redistribution in the sample which makes the sample 
edge become negatively charged. In the vortex state there is a competition between this Meissner 
charge and the vortex charge which may change the polarity of the charge at the sample edge with 
increasing magnetic field. It is shown analytically that in spite of the charge redistribution the 
mesoscopic sample as a whole remains electrically neutral. 
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C/3 I. INTRODUCTION 



Recently it was predicted that the core of an Abrikosov vortex in bulk type-II superconductors is charge d&l. This 
effect occurs because of the difference of the chemical potential in the superconducting versus normal state. Such 
a change in chemical potential results in a redistribution of the electrons in the regionjiear the vortex core and 
£^ ■ culminates in a charging of the vortex core when the superconductor is in the mixed stateao. 

In the present paper we investigate this phenomenon in mesoscopic superconductors. A mesoscopic sample has 
a typical size which is comparable to the coherence length (£) and the magnetic field penetration length- (A). The 
behaviour of such structures in an external magnetic field (H ) is strongly influenced by the sample shapeQU and may 
y—i , lead to various superconducting states and different phase transitions between them. Jumps in magnetization were 
observed when the applied magnetic field or temperature (T) are variedS. 

A number oteaxlier works studied the geomet ry , d ependent magnetic respo nse of mesoscopic superconductors: i) disk 
shape samples0i!3; ii) infinitely long cylindersEJO; iii) ring-like structurestL30 and more complicated geometries^ 8 ! 



Theoretical studies of mesoscopic superconductors are based on the phenomenological Ginzburg-Landau (GL) theor 
I i which successfully describes mesoscopic samples in a wide H-T region. In particular, it has been shown that in 
mesoscopic samples (disks or cylinders) surrounded by a vacuum or an insulator two kinds of superconducting states 
can exist. Firstly, there is a circular symmetric state with a fixed value of angular momentum, called giant vortex. 
The observed magnetization jumps correspond to first order phase transitions between giant vortices-.with different 
angular momentumEH. Secondly, in samples with a sufficiently large radius multi-vortex structures^ can nucleate 
which are the analogue of the Abrikosov flux line lattice in a bulk superconductor. These states can be represented 
as a mixture of giant vortex states with different angular momentum. For multi-vortex states it is also possible to 
introduce an effective total angular momentum, which is nothing else then the number of vortices in the disk, i.e. the 
vorticity. With cha*wieg the magnetic field there is a second order phase transition between the multi-vortex and the 
Q , giant vortex state@li3'EI 

It is expected, that the charge distribution in such mesoscopic samples may be appreciably altered due to the pres- 
ence of a boundary. Furthermore, screening currents near the boundary of the sample will also lead to a redistribution 
of charge and consequently even in the Meissner state there will be a non-uniform distribution in the sample. In a 
certain sense, this case can be viewed as a vortex turned inside out, i.e. with its core at infinity. In the presence of 
vortices there will be an interplay between the Meissner charge and the previously studiecU~0 vortex charge. 

The present paper is organized as follows. In Sec. II we give the necessary theoretical formalism on which our 
numerical results arc based. In Sec. Ill we investigate the charge distribution in both the Meissner state and the giant 
vortex states. Then the charge distribution in the multi- vortex state is discussed for thin disks (Sec. IV). Our results 
are summarized in Sec. V. In the Appendix we present the proof of electrical neutrality in a mesoscopic superconductor 
of general shape. 
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II. THEORETICAL APPROACH 



We consider a mesoscopic superconducting sample of circular symmetry with radius R and thickness d surrounded by 
an insulating medium. The external magnetic field H = (0, 0, H) is uniform and directed normal to the superconductor 
plane. The starting point of our analysis is that the rotating motion of Cooper pair* around the vortex core leads to 
a spatial redistribution of charge carriers, which generate the electrostatic potentialTH 

^(T > ) = ^ (|V'(^)| 2 -l), (1) 

where (~f*) is the dimensionless superconducting order parameter normalized so that \if> ("r*)) 2 is measured in units 
of the Cooper pair density in a bulk superconductor. The amplitude ipo is different in different approaches. We use 
fo xi [ckJ /2e as proposed in Ref.B. A three-times smaller value tpo — \a\ /6e was used in Ref.EJ. From the theory in 
Rcf.cl it approximately follows <po = ( | ck| /27r 2 e) (dT c /dln.€F ), where dT c /dlneF ~ In (HuJo/kBTc) ~ 1 to f 0. Thus all 
approaches yield an electrostatic potential of a similar magnitude. r— ■ 
The distribution of the corresponding charge density q (T*) is obtained from the Poisson equationtll 

Anq {?) = ip (!*) . (2) 



The Cooper pair density \ip (~r*)\ 2 is determined from a solution of the system of two coupled non-linear GL equations 
for the superconducting order parameter, ip{~i*), and the magnetic field (or vector potential A(~r*)) 

(-^-^)% = ^»-VM 2 , (3) 



K 2 V X V X ~A 



(4) 



where the density of the superconducting current j is given by 



(•5) 



Here T 5, = (~p,z) is the three-dimensional position in space. Due to the circular symmetry of the sample we use 
cylindrical coordinates: p is the radial distance from the disk center, 9 is the azimuthal angle and the z-axis is taken 
perpendicular to the disk plane, where the disk lies between z = —d/2 and z = d/2. For d — > oo we obtain the 
cylinder geometry. All distances are measured in units of the coherence length £ = H/ y2r 
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| a | (m* = 2m is the mass 

of the Cooper pair), the vector potential in ch/2e£, the magnetic field in H C 2 — ch/2eS l 2 = k^/2H c , where H c is the 
thcrmodynamical critical field, the superconducting current in jo = cH c /2tt^, and n = A/^ is the GL parameter. 

Eqs. (|H|) has to be,supplemented by boundary conditions for ip(!*) and ~A(T f ). For the superconducting condensate 
it can be written ast2l 



(-iV - it) ip 



0, 



(6) 



where ~ft is the unit vector normal to the sample surface. The boundary condition for the vector potential has to be 
taken far away from the superconductor where the magnetic field becomes equal to the external applied field H 



A 



1 



-Hp~t g , 



(7) 



where itg denotes the azimuthal direction. 

The free energy of the superconducting state, measured in F = H^V/8tt units, is determined by the expression 



F 



2 
V 



dV 



> 4 



(8) 



with the magnetic field 



~t = V x ~A{T > ). 
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Here we will consider three different geometries for the superconductor: I) an infinitely long cylinder; II) a thin disk 
with finite thickness d < A,£; and III) an infinitely thin disk, i.e. d— > 0. In all three cases the superconducting order 
parameter does not depend on z. In case I it obviously follows from the sample geometry. In cases II and III it was 
foundtj that the dependence of i/j(~f*) on z is very slow. This allows us to average ip(~J*) over the sample thickness and 
to solve the problem for the two-dimensional problem for the order parameter if) (p, 9). However, the magnetic field 
in case II has a z-dependence, which is responsible for the demagnetization effect. For both cases I and II we soli^e 
the problem numerically by the method proposed in Ref.B. For thin mesoscopic disks we use the results of Refs. &E2I 
which allowed us to solve the problem semi-analytically. 



III. CHARGE IN THE MEISSNER AND THE GIANT VORTEX STATES 



Firstly, we consider the situation with a fixed value of the vorticity L. The giant vortex state has cylindrical 
symmetry and consequently the order parameter can be written as ip (~f?) = f (p) exp (iL9). For a thin disk (case III)) 
the order parameter isE3 

T \ V 2 



VG>,0)=(-A^) f L (p)exp(iL9), (9) 



where 

/i(p)=( r ^l expl-^lMl-i/.L + l,^), (10) 



Hp 2 \ L/ \ ( B(?\ / r , , Up 2 



h = J Pdp ft (p) , h = J Pdp fl (p) , 



(11) 



k = H (l + 2i/) - 1. 



(12) 



Here M (a, b, y) is the Kummer functioned and the value of v is determined by the non-linear equation, which results 
from the boundary condition (H) 



L ) M ( —v, L + l,— , 

2 J V 2 J L + l 



M I _i/ + l,i + 2 > - 



0. 



(13) 



Here $ = HR 2 is the magnetic flux through the disk in the absence of any flux expulsion. Using (p|-|l"2|) we can derive 
explicitly the charge distribution 



h fH 



q { P )=±Kf[-) p 



„2(£-l) 



exp 



Hp 2 



{ L-M- { -v, L + l, - (2L + 1) M (-v, L+l,^- 



M —v, L + l 



2 r 



Hp 2 



2v ( Hp 2 
-M [ - u + l,L + 2,-g- 



(Hp 2 



V 2 

xM (-v+l,L + 2, 
2v(v+l) 



M- | —v. L + 1. 
tfp 2 



L + l 
LTp 2 



2 / L + l 

2^ 2 



Hp 2 



-v + l,L + 2, % 

(L+l) 2 V '2 



1/ i - l /,L + l,^Mf-H2,L + 3 



(14) 



(L + l)(L + 2) 

where the charge density is measured in units of qo = ti 2 / (167r?7i*e£ 4 ) = (as/327r£ 4 ) e (as — h 2 /me 2 is the Bohr 
radius). For Al we have £ = 250 nm □ which results into qo/e ~ 1.3 • 10~ 13 nm -3 . 
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A. MEISSNER STATE 



First, we study the Meissner state, i.e. L — 0. The radial dependences of the Cooper pair density (or the corre- 
sponding distribution of the potential <p (p) (see Eq. (0)), the charge density q(p) and the density of superconducting 
current j (p) are shown in Figs. Ill(a-c) for the cases of cylinder, finite disk and very thin disk, respectively, for different 
values of the applied field. Fig. fl[(d) shows these dependencies for a cylinder at the field H = 0.52i/ C 2 but for different 
k values. Due to the finite radial size of the samples all distributions are inhomogeneous along the radius of the 
sample. The Cooper pair density is maximum at the center and decays towards the sample edge. As a result, in 
the center of the sample there is a region of positive charge while near the edge a negative "screening" charge q scr is 
created. To avoid confusion let us note that for simplicity we write "positive charge" instead of "charge of the same 
sign as is the sign of the dominant charge carriers" . In the cylinder the Cooper pair density decreases with increasing 
field, while both the screening superconducting current and the charge polarization monotoneously increase. The 
behavior for the disks is more complicated. For small fields the picture is very similar to the one of the cylinder. 
But for fields where the Meissner state becomes metastable (i.e., H/H C 2 > 0.32 for the parameters used in Fig. |f|(b)) 
the screening charge becomes maximal. The positive charge is pulled further towards the center of the disk and the 
screening charge region expands. Notice that now the maximum of the positive charge decreases with field and also 
at the surface its absolute value decreases. For a cylinder, at fixed field but with increasing k (Fig. |l|(d)) its charge 
distribution behaves similar as for fixed n and increasing magnetic field (see Fig. 0(a)). 

Notice that even for L = when no vortex is present inside the superconductor, there is still a non-uniform charge 
distribution. This charge redistribution can be characterised by two quantities: i) the distance p* which separates the 
positive and negative charge regions, and ii) the total screening charge Q-^ scr which is defined by the integral over 
the region VL occupied by the negative charge: 

Q-, SC r = 2nd J pq scr (p)dp. (15) 

(p-) 

The dependences of p* (H) are shown in Fig. ||(a-c) for the cylinder, the finite disk and the thin disk, respectively. In 
Fig. |(c) the open squares refer to the magnetic field region where the L = state is metastable (the crossed circles will 
be explained below). Notice that p* decreases with increasing field and this decrease is more pronounced for thinner 
disks. The magnetic field dependences of the absolute value \Q—\ of the screening charge is shown in Fig. ||(a-c) for 
the same geometries. The screening charge Q SC r increases with magnetic field, but for the disk cases a local maximum 
is reached. This local maximum is reached for fields where p* starts to decrease more strongly. Taking £ = 250 nm 
and as = 0.05 nm we obtain that for the characteristic value of screening charge |Q_ ; scr \ ~ 3 • 2n^ 2 q d ~ 8 • 10 _6 e 
at the field H = Q.3H C 2 for the Al disk with d/£ = 0.2 (the case II). For strongly type-II superconductors which have 
a very small coherence length £ this induced charge can be made orders of magnitude larger. 

By direct integration of q (p) over the sample surface one can convince ourselves that the total charge per unit of 
sample length Q = § d9 f Q pq (p) dp = 0, i.e. there is charge neutrality over the whole sample. But this fact can be 
easy generalized analytically for states with any vorticity. The proof is given in the Appendix. 



B. GIANT VORTEX STATE 

The same dependences as in Figs. fj](a-d) are shown in Figs. |](a-d) for samples in the L = 1 vortex state. In this 
state \4> (p)\ 2 — in the center of the vortex core located in the centre of the sample. Notice, that for a cylinder the 
charge distribution almost does not change with magnetic field. The reason is that the external field affects the Cooper 
pair density only near the sample edge region. This is different for the disk geometry, where large demagnetization 
effects strongly influences the penetration of the magnetic field in the disk which changes the vortex structure. The 
vortex core is negatively charged and at small fields the positive charge outside the vortex extends up to the border 
of the sample. The dependences of the position p* (H) and the size of the charge pile- up \Q- \ (H) for the L = 1 state 
are shown in Figs. ||(a-c) and |](a-c). 

With increasing external magnetic field the screennig current at the sample surface increases and makes the region 
near the border of the sample negatively charged. In this case there exist two p* , where q(p*) = 0. The core of 
the vortex is negatively charged with total charge \Q~. „|. Around this core there is a ring of positive charge which 
compensates the charge of the vortex core and the surface charge. Near the surface a ring of negative screening charge 
exists with total charge \Q— S cr\- The size of the latter increases with increasing magnetic field. For disks, \Q~, S cr\ 
reaches a local maximum after which it decreases in the large magnetic field region; this is the region where the L = 1 
state is unstable (see circles with crosses in Fig. ||(c)). 
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Next we investigated the charge distribution for the vortex states with L > 1 and we limited ourselves to the thin 
disk case. From Eq. ([l4| ) one finds immediately that q — in the center of the disk when L > 2. Consequently, the 
charge distribution in the vortex core has a ring shape. To illustrate this we show in Figs. ||(a,b) the same dependences 
as in Fig. ||(c) but now for L = 2 and 3, respectively. With the exception of the core region the charge distribution 
for the giant vortex states is qualitatively similar to the case L = 1, and the charge on the sample surface changes 
sign with increasing external magnetic field. Notice also that the number of areas where the charge changes its sign 
does not increase with L: it equals 2 for small fields and increases to 3 for higher magnetic fields. 

IV. VORTEX CHARGE IN THE MULTI- VORTEX STATE 

For sufficiently large radial size of the superconductor the giant vortex state can-break up into multi-vorticeJHll. 
To explain the physics we limit ourselves to the case of a thin disk. It was shownEjO that the order parameter of the 
multi-vortex state in general can be viewed as a superposition of giant vortex states with different Lj 



V>(7?)= E C L J L] (p)exp( l L J 9). 



Lj=0 



(16) 



where L is now the value of the effective total angular momentum which equals the number of vortices in the disk. 
For disks with not so large radius (3 5 £) the. order parameter of the multi- vortex state is the superposition of only 
two states and is described by the expressionEa 



ip (j?) = C L J Ll (p) exp {iL x 9) + C L2 f La (p) exp (iL 2 9) , 



(17) 



where 



Cr 



—Al 1 Al 2 Bl 1 + 2Al 2 Al 1> l 2 Bl 2 



A L ,A L .-4A 2 r 



Li,L 2 



—Al 2 A Li B L2 + 2A Ll A LltL2 B Ll 
A Ll A L2 - AAl i L2 



1/2 



1/2 



(18) 



A, 



2nd 
~V~ 



Pdp Shi 0°) ' 



2nd 
~V~ 



IS 

J pdp si ( P )Sl 2 (p), b Li 



2nd 
~V~ 



P d P fit (P) ) 



(19) 



and Shi (p) and A^ 4 are determined by Eqs. (|10|) and (|12[), respectively. The charge density distribution is then given 
by the expression 



q(p,0) =2C 2 Ll 



S'h\ (P) + f Ll (p) S'L (P) + -hi (P) S'hr (P) 

p 

St (p) + h> (p) f'L (p) + (p) S'h 2 (p) 



(20) 



-2C Ll C L2 cos [{L, - L 2 ) 9} [f^ (p) f L2 (p) + 2f' Li (p) fh 2 ( P ) + f Ll (p) S'L (p) 



+-(S Li (p)Sl 2 ( P ) + f Ll (p)S'l 2 (p)) 



{Li - L 2 y 



S^ (p)Sh 2 (p) 



where the prime denotes the derivative with respect to p. The explicit expression is rather lengthy and is therefore 
not given here. ___ 

Earlier analyses have showriBESEJ that there exist two kinds of multi-vortex states: i) stable configurations which 
correspond to a minimum of the free energy; and ii) states which correspond to saddle points of the free energy. 
The latter ones correspond to the energy barrier states between states with different vorticity L and describe the 
penetration of flux into the disk. Due to the transitions between the different L states with increasing (or decreasing) 
external field some giant vortex states are never realised (for example, such states for L = and 1 correspond to the 
crossed circles in Figs. ||(c) and ||(c)). 
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As an example, we consider a thin disk with R/£ — 4.0. The charge density q(x,y) distribution over the disk for 
the different kinds of multi- vortex states are shown in Figs. |^ and]?]. In Fig. ^ this distribution is given for the saddle 
point state (1:2) at the field H = 0.32if C 2 and the contour plot of the distribution of the Cooper pair density |^| for 
this state is shown at the bottom of the figure. The dark regions on the contour plot correspond to low Cooper 
pair density. The same distributions for the stable multi-vortex state (0:4) at the field H = 0.75_ff C 2 are shown in 
Fig. [?]. One can see region s of negative charge located at the vortex core and positive charge near the edge of the 
sample. Using Eqs. ( |P7H1S| ) it is easv to prove that also in the multi-vortex state the disk is electrically neutral as a 
whole (see also the Appendix). 



V. CONCLUSIONS 



We studied theoretically the redistribution of electrical charge in circular mesoscopic superconducting samples with 
different shape, i.e. disks and cylinders. The theory applies for intermediate temperatures where the Ginzburg-Landau 
theory still gives reasonable results while the share of the superconducting electrons is already of order unity so that 
the screening by normal particles may be neglected. Previously, it was predicted that the vortex core in bulk type-II 
samples is negatively charged. Here we found that even in the Meissner state with no vortices inside the sample 
there exists a non-uniform charge distribution. Due to the finite radial size a region near the sample edge becomes 
negatively charged while the interior of the sample has a corresponding positive charge. This charge redistribution 
is a consequence of the screening currents near the sample edge which makes it behave like a vortex which is turned 
inside out. When vortices are inside the sample there is a superposition of the vortex charge and this Meissner charge. 
Because of this interplay between vortex charge, which is positive near the sample surface, and the Meissner charge, 
which is negative at the sample surface, the charge at the sample edge changes sign as a function of the applied 
magnetic field. We also proved analytically that the there is only a redistribution of charge and that the total sample 
charge is neutral as long as the boundary condition (0) is satisfied. 
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VII. APPENDIX: PROOF OF ELECTRICAL NEUTRALITY IN A MESOSCOPIC SAMPLE 



Using Eqs. (|l|j^) and Gauss theorem the total charge Q can be expressed as 



(Al) 



2e 

H 

2e 



ip* ( r») VV (~?) + "0 {~?) vV* ("?*) 



dS 



ip* ( r>) - il^j iP{~r>) + i}) ( r>) ( V* + iTj i>* ( ~?) 



dS. 



From the boundary condition (^|) and its complex conjugate it follows immediately that Q = 0. Notice, that this 
result is very general, it is valid for any vortex configuration (the giant vortex states and the multi-vortex ones) and 
arbitrary shape of the superconducting sample as long as the boundary condition (||) is satisfied. 
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FIG. 1. The radial dependences of the Cooper pair density |-;/;(p)| 2 , the charge density q (p) and the supercurrent density 
j (p) in the Meissner state for (a) an infinite long cylinder with k = 1.0, (b) finite thickness disk with d/£ = 0.2, k = 1.0, 
(c) very thin disk for different magnetic fields, and (d) for an infinite long cylinder with different Ginzburg-Landau parameter 
K at the magnetic field H — 0.52//^. All samples have the same radius R/£ = 4.0. 

FIG. 2. The position of the boundary between the regions of negative and positive charges inside the sample as a function of 
the external magnetic field for: (a) the infinite long cylinder with k = 1.0, (b) the finite thickness disk with d/£ = 0.2, k = 1.0, 
and (c) the thin disk. The p* (H) dependences for the Meissner state are shown by the solid line, and for the vortex state - by 
the dashed lines. The dash-dotted curves represent the position of zero current. 

FIG. 3. The absolute value of the negative charge inside the sample as a function of the external magnetic field for: (a) the 
infinite long cylinder with n — 1.0, (b) the finite thickness disk with d/£ = 0.2, n — 1.0, and (c) the thin disk. The dependences 
for the Meissner state are shown by the solid line, and for the vortex state - by the dashed lines. 

FIG. 4. The same as in Fig. [j] but for the single vortex state (with L — 1). 

FIG. 5. The radial dependences of the Cooper pair density |-!/>(p)| 2 , the charge density q (p) and the supercurrent density 
j (p) in the giant vortex state with (a) L — 2 and (b) L = 3 for a thin disk with 7?/£ = 4.0 and for different applied magnetic 
fields. 
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FIG. 6. The charge density distribution for the (1:2) saddle-point state at the field H = 0.32f/ c2 for the thin disk with 
Jl/£ = 4.0. The bottom contour plot shows the distribution of the Cooper pair density. 

FIG. 7. The same as in Fig. |but for the (0:4) multivortex state in the disk with R/£ = 4.0 at the field H = 0.75H c2 - 
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